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Abstract 

It is shown that the condition for chemical freeze-out, average energy per hadron 
~ 1 GeV, selects the softest point of the equation of state, namely the point where 
the pressure divided by the energy density, p[e)/e has a minimum. The sensitivity 
to the equation of state used is discussed. The previously proposed mixed phase 
model, which is consistent with lattice QCD data naturally leads to the chemical 
freeze-out condition. 
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Over the last few years the question of chemical equilibrium in heavy ion collisions 
has attracted much attention Assuming thermal and chemical equilibrium within 
a statistical model, it has now been shown that it is indeed possible to describe the 
hadronic abundances produced using beam energies ranging from 1 to 200 AGeV. The 
observation was made that the parameters of the chemical freeze-out curve obtained at 
CERN/SPS, BNL/AGS and GSI/SIS all lie on a unique freeze-out curve in the T — hb 
plane. Recently, a surprisingly simple interpretation of this curve has been proposed: the 
hadronic composition at the final state is determined solely by an energy per hadron of 
approximately 1 GeV per hadron in the rest frame of the system under consideration 
P, ^. In this letter we propose a dynamical interpretation of the chemical freeze-out 
curve and show that it is intimately related to the softest point of equation of state 
defined by the minimum of the ratio -(e) as a function of e f^. Our considerations are 
essentially based on the recently proposed |^ mixed phase model which is consistent 
with the available QCD lattice data 0. The underlying assumption of the Mixed Phase 
(MP) model is that unbound quarks and gluons may coexist with hadrons forming an 
homogeneous quark/gluon-hadron phase ||. Since the mean distance between hadrons 
and quarks/gluons in this mixed phase may be of same order as that between hadrons, 
their interaction with unbound quarks/gluons plays an important role defining the order 
of the phase transition. 

Within the MP model ^ the effective Hamiltonian is written in the quasi particle 
approximation with the density-dependent mean-field interaction. Under quite general 
requirements of confinement for color charges, the mean-field potential of quarks and 
gluons is approximated by the following form: 

U,{p) = U,{p) = ^ (1) 

with the total density of quarks and gluons 

P = Pq + Pg + Y. ^JPj 
j 

where pg and pg are the densities of unbound quarks and gluons outside of hadrons, 
while pj is the density and nj is the number of valence quarks inside the hadron of 
type j. The presence of the total density p in (|l|) corresponds to the inclusion of the 
interaction between all components of the mixed phase. The approximation (0) recovers 
two important limiting cases of the QCD interaction, namely, if p the interaction 
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potential goes to infinity, i.e. an infinite energy should be spent to create an isolated 
quark or gluon which ensures the confinement of color objects and, in the other extreme 
case of high energy density corresponding to p oo we obtain the asymptotic freedom 
regime. 

The use of a density-dependent potential for quarks and a hadronic potential 
described by a modified non-linear mean-field model requires certain constraints, re- 
lated to thermodynamic consistency, to be fulfilled 0, For the chosen form of the 
Hamiltonian these conditions require that Ug{p) and Ug{p) should be independent of the 
temperature. From these conditions one also obtains an expression for the form of the 
quark-hadron potential p. 

A detailed study of the pure gluonic SU{3) case with a first order phase transition 
allows one to fix the values of the parameters as 7 = 0.62 and A^/^^^+^^ = 250 Mel/. 
These values are then generalized to the the SU (3) system including quarks. For the case 
of quarks of two light flavors at zero baryon density, ub = 0, the MP model is consistent 
with the results from lattice QCD with a deconfinement temperature Tdec =153 MeV 
and the crossover type of the deconfinement phase transition. The model can be extended 
to baryon- rich systems in a parameter-free way [Q. 

A particular consequence of the MP model is that for hb = the 'softest point' of 
the equation of state, as defined in Q], is located at comparatively low values of the 
energy density: esp ~ 0.45 GeV / fm? . This value of e is close to the energy density 
inside a nucleon and, thus, reaching this value signals us that we are dealing with a single 
'big' hadron consisting of deconfined matter. For baryonic matter the softest point is 
gradually washed out at hb > 0.4 hq. As shown in [^,0, this behavior differs drastically 
from both the interacting hadron gas model which has no soft point and the two-phase 
approach, based on the bag model, having a first order phase transition by construction 
and the softest point at esp > 1 GeV/ frn? independent of hb @]. These differences 
should manifest themselves in the expansion dynamics. 

In Fig.l we show trajectories of the evolution of central Au-|-Au collisions in the T — ^b 
plane together with the freeze-out parameters obtained from hadronic abundances. The 
initial state was estimated using a transport model starting from a cylinder in the center- 
of-mass frame with radius R = 4 fm and length L = 2R/'yc.m. as described in 0, 0. The 
subsequent isoentropic expansion was calculated using a scaled hydrodynamical model 
with the MP equation of state. As seen from the figure, the turning points of these 
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trajectories correlate nicely with the extracted freeze-out parameters, as was noted in 0, 
as well as with the smooth curve corresponding to a fixed energy per hadron in the 
hadronic gas model [Q]. 

The observed correlation is further elucidated in Fig. 2. The quantity p/e is closely 
related to the square of the velocity of sound and characterizes the expansion speedQ, so 
the system lives for the longest time around the softest point which allows it to reach 
chemical equilibrium for the strongly interacting components. It is also seen that the 
position of the softest point correlates with the average energy per hadron being about 
1 GeV in all nuclear cases and even for pp collision. One should note that the quantity 
^1 Phad, where e is the total energy density, coincides with {Ehad) / {Nhad) considered in 
only in the case when there are no unbound quarks/gluons in the system. In the MP 
model, all components are interacting with each other and therefore the quantity (Ehad) 
is not defined. The admixture of unbound quarks at the softest point esp amounts to 
about 13% and 8% at beam energies Eiab = 150 and 10 AGeV, correspondingly. 

The MP equation of state plays a decisive role for the regularity considered here, 
describing both the order of the phase transition and the deconfinement temperature. 
The two-phase (bag) model exhibits a first order phase transition with T^ec = 160 MeV 
and has a spatially separated Gibbs mixed phase but the corresponding trajectories in 
the T — fiB plane are quite different from those in the MP model as shown in 0. The 
exit point from the Gibbs mixed phase at Eiab = 150 AGeV is close to the corresponding 
freeze-out point in Fig.l. However the large differences noted above in esp and in its 
dependence on n^, mainly caused by the different type of the predicted phase transition, 
does not lead to the observed correlation with the softest point position in the whole 
energy range considered. The interacting hadron gas model has no softest point effect as 
was demonstrated in 0, |^. This fact is seen also from Fig. 2 where at Eiab = 2 AGeV 
the quark admixture is practically degenerated (~ 1%) and instead of a minimum there 
occurs a monotonic fall-off specific for hadronic models with a small irregularity in p/e 
near the point e/phad = 1 GeV f\. 

It is noteworthy that similarly to the results presented in Fig. 2 the softest point of 
the equation of state correlates with an average energy per quark, e/p ~ 350 MeV 

*In simple hydrodynamic models, for example, the transverse expansion of a cylindrical source, the 

evolution is governed by the pressure-to-enthalpy ratio, p/{p + e) [lO[ |. 

tNote that at the SIS energies the chemical freeze-out point practically coincides with the thermal 
freeze-out 
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which is close to the constituent quark mass. So, at higher values of e/ p we are deal- 
ing with a strongly-interacting mixture of highly-excited hadrons and unbound massive 
quarks/gluons forming (in accordance with Landau's idea |]12|) an 'amorphous' fluid suit- 
able for hydrodynamic treatment. Below the soft point the interaction deceases, the 
relative fraction of unbound quarks/gluons decreases, higher hadronic resonances decay 
into baryons and light mesons and thereby the value of e/p goes down. 

In summary, the unified description of the chemical freeze-out parameters found in 
is naturally related to the fact that the proposed condition {Ehad) / {Nhad) ~ 1 GeV 
selects the softest point of the equation of state where the strongly interacting system 
stays for a long time. Such a clear correlation is observed for the equation of state of the 
mixed phase model but not in purely hadronic nor in two-phase models. In this respect 
the success of the MP model in the dynamical interpretation of the freeze-out regularity 
may be considered as an argument in favor of a crossover type of the deconfinement phase 
transition in 5*^7(3) system with massive quarks. 
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Figure captions 



Fig.l. The compiled chemical freeze-out parameters (borrowed from |^) obtained 
from the observed hadronic abundances and dynamical trajectories calculated for central 
Au + Au collisions at different beam energies Eiab with the mixed phase equation of state. 
The smooth dashed curve is calculated in the hadronic gas model for (Ehad) / {Nhad) = 
IGeV i. 

Fig. 2. The ratio of pressure to energy density, p/e, versus the average energy per 
hadron, e/phad, for evolution of different systems. The upper curve corresponds to pp 
collisions at = 40 GeV with isoentropic expansion from a sphere with R = 1 fm. 
Other cases are calculated for central Au + Au collisions at the given beam energy under 
the same conditions as in Fig.l. 
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